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Magnetic fields in galaxies are produced via the amplification of seed magnetic 
fields of unknown nature. The seed fields, which might exist in their initial 
form in the intergalactic medium, were never detected. We report a lower 
bound B > 3 x 10^^ gauss on the strength of intergalactic magnetic fields, 
which stems from the nonobservation of GeV gamma-ray emission from elec- 
tromagnetic cascade initiated by tera-electron volt gamma-ray in intergalactic 
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medium. The bound improves as ' if magnetic field correlation length, A s, 
is much smaller than a megaparsec. This lower bound constrains models for 
the origin of cosmic magnetic fields. 

The problem of the origin of 1- to 10-yuG magnetic fields in galaxies and galaxy clusters is 
one of the long-standing problems of astrophysics and cosmology [see {1-4) for reviews]. It is 
assumed that the observed magnetic fields result from the amplification of much weaker seed 
fields. However, the nature of the initial weak seed fields is largely unknown. There are two 
broad classes of models for the seed fields: astrophysical models, which assume that the seed 
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fields are generated by motions of the plasma in (proto)galaxies, and cosmological models in 
which the seed fields are produced in the early universe {1-4). 

Extremely weak unamplified extragalactic magnetic fields (EGMFs) have escaped detection 
up to now. Measurements of the Faraday rotation in the polarized radio emission from distant 
quasars (7, 5, 6) and/or distortions of the spectrum and polarization properties in the cosmic 
microwave background (CMB) radiation (7-18) imply upper limits on EGMF strengths at the 
level of ~ 10~^ G. Numerical modeling of magnetic field formation in galaxy clusters implies 
a theoretical upper bound of the order of ~ 10~^^ G on EGMF strength (19, 20). Bounds on 
the EGMF strength depend on the field correlation length A^, which is also unknown. A lower 
limit on As is set by the requirement that the resistive magnetic diffusion time scale has to be 
larger than the age of the Universe (2), whereas an upper limit is set only by the size of the 
visible part of the Universe, Rh- 

Here we report a lower bound for the EGMF strength, derived from the data of Fermi and 
High Energy Stereoscopic System (HESS) gamma-ray telescopes. Similarly to the existing 
upper bounds, the lower bound depends on the unknown EGMF configuration parameters, such 
as the typical correlation length and spectrum. 

Gamma rays with energies above ~ 1 TeV cannot propagate over cosmological distances 
because of absorption resulting from interactions with diffuse extragalactic background light 
(EEL) (21-24). The mean free path of gamma rays of energy E^^ through EBL is ~ 
80ft; {E^JIO TeV)"' Mpc, where ft; ~ 1 is a numerical factor that accounts for uncertainties 
of the measurements and modeling of the EBL (25). Interactions of multi-TeV gamma rays 
with the EBL lead to the deposition of electron-positron pairs in the intergalactic space. These 
e+e~ pairs emit secondary cascade gamma rays via Inverse Compton (IC) scattering of CMB 
photons. Typical energies for the IC photons emitted by electrons of energy ~ Ejg/2 are 
E^ = {4:/3)ecMB{Ee/meC^f ~ SS[E^J10TeVf GeV, where ecmb = 6 x 10"^ eV is the 
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typical energy of CMB photons and Ee and rUe are the energy and mass, respectively, of an 
electron. Pairs lose energy on IC scattering on the distance scale Dg ~ 10^^(i?e/10 TeV)^^ cm, 
which is much smaller than the gamma ray mean free path D^. Power removed from the primary 
gamma-ray beam is transferred to the cascade gamma-ray emission. 

If magnetic fields, which deviate electron and positron trajectories, are negligibly small, 
the IC emission from the electromagnetic cascade contributes to the primary point gamma-ray 
source flux (26-29, 54). Otherwise, if magnetic fields along the path of development of the 
cascade are strong enough to deviate the trajectories of the pairs, the cascade emission appears 
as extended emission around the initial point source (25, 30-33). 

The deflection angle 5 depends on the correlation length of the magnetic field. If As » De, 
the motion of electrons or positrons can be approximated by the motion in a homogeneous mag- 
netic field. In this case 5 ~ D^/Rl ~ 3 x 10"^ [B/IQ-^^ G] [EjlQ TeV]"^ is a ratio of to 
the Larmor radius i^L. If As <S D^, electron deflections are describable by diffusion in angle, so 
that the deflection angle is 5 = ^fDj^jRi^ ~ SxlO'^ [Ee/10 TeV]"^/^ [5/10-^^ G] [As/l kpc]^/^ 
The size of the extended cascade source is estimated as Gext — 5/t, where r — D/D^is the 
optical depth for gamma rays from a source at a distance D with respect to absorption on 
EBL (25). 

Because lower energy electrons are deviated by larger angles, the size of the extended cas- 
cade source ©ext is larger at low energies. The energy of cascade photons, Ej^jnin-, below which 
the extended source size becomes larger than the point spread function (PSF) of a telescope, de- 
pends on the EGMF strength and correlation length. In the case of the Fermi telescope, the PSF 
depends on the photon energy, decreasing as ©psf — 2° [E^/1 GeV]~°'^ (95% of the signal) 
below EjC^l GeV and improving from ~ 2° at 1 GeV to ©psf ^ 0.2° at E^ ~ 10 GeV (34). 
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Taking the photon energy E. 



10 GeV as a reference, one finds that 6ext > ©psf if 



B> B 



PSF — 



6 X [E^,min/10 GeV] G, 

8 X lO-i^r [E^,min/10 GeVf" [Xb/1 kpc] 



' G, Xb < De 



(1) 



To constrain the presence of a cascade contribution in the spectra of distant TeV blazars, 
we have analyzed Fermi/LAT (Large Array Telescope) data on the blazars lES 1101-232, 
lES 0229+200, lES 0347-121 and H 2356-309, obtained during the Fermi's first year of 
operation. These sources were selected on the basis of their high redshifts (z = 0.14 for 
lES 0229-1-200, z = 0.165 for H 2356-309, z = 0.186 for lES 1101-232 and z = 0.188 
for lES 0347-121) and hard TeV band spectra (22, 23, 35-38). 

Minimal possible cascade signal from lES 1 101-232 and H 2356-309 are below the Fermi 
upper bounds [see supporting online material (SOM) for details] . The nondetection of lES 0229-1-200 
and lES 0347-121 with Fermi imposes restrictions on the contribution to the flux from the low- 
energy tail of the cascade (40). In order to determine these restrictions, we have developed a 
numerical model of electromagnetic cascade in intergalactic space, initiated by the absorption 
of TeV gamma rays on EBL photons. The model solves the one-dimensional kinetic equations 
for concentrations of gamma rays and of e"'"e^ pairs, taking into account the injection of pairs by 
the absorbed gamma rays (41), the injection of secondary cascade gamma rays via IC emission 
by the e~^e~ pairs and the cooling of the e~^e~ pairs via IC scattering energy losses (42). The 
cosmological photon backgrounds with which our gamma rays and e+e" pairs interact include 
CMB and EBL taken from (24). 

The initial gamma-ray spectra at the source, shown by dashed thin curves in Fig. [H were 
chosen in the form of a cutoff powerlaw, dN^/dE ~ E~^ exp{—E / E^ut)- The cascade emis- 
sion power is equal to the fraction of the power of the primary gamma-ray beam absorbed on 
the way from the source to Earth. Because almost 100% of the power initially injected at the 
energies above TeV is absorbed, the luminosity of cascade emission is roughly equal to the in- 
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Fig. 1 : A comparison of models of cascade emission from TeV blazars (thick solid black curves) 
with Fermi upper limits (grey curves) and HESS data (grey data points). Thin dashed curves 

show the primary (unabsorbed) source spectra. Dotted curves show the spectra of electromag- 
netic cascade initiated by pair production on EBL. Vertical lines with arrows show the energies 
below which the cascade emission should be suppressed. 
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tegral primary source luminosity in the multi-TeV energy band. Primary source spectra shown 
in Fig. [TJcorrespond to the minimal value of Ecut compatible with HESS data, to minimize total 
flux of the cascade contribution (see details in the SOM). 

Assumption of zero magnetic field along the cascade development path is in contradiction 
with Fermi upper bounds on the source fluxes (Fig. [T]). Although the model spectra deviate from 
simple powerlaws, the deviations are small, meaning that Fermi bounds on the powerlaw-type 
spectra could be applied. 

The cascade emission has to be suppressed below an energy E^^min (marked by a vertical 
line in the three graphs of Fig. [T]) at which the Fermi upper bound becomes higher than the 
model cascade flux. Suppression of the cascade flux at low energies could be achieved if tra- 
jectories of low-energy e+e~ pairs are deviated by magnetic fields. The cascade emission does 
not contribute to the point source flux below the energy E-y^rnin if EGMF is stronger than -BpsF 
given by Eq. 1. The values of -BpsF corresponding to E^^min found for each source are given in 
table SI. 

The best bound B > BpsF, imposed by Fermi limits on the flux from lES 0229-1-200 (black 
hatched region in Fig. |2l), suffers from a number of uncertainties and, therefore, should be 
considered as an order-of-magnitude estimate. 

Suppression of the cascade contribution to the point source flux below the energy E^^rnin 
results in a deviation of the model source flux from the powerlaw at the energies E < Ey^min- 
This means that the Fermi upper bounds on the powerlaw-type spectra shown in Fig. [T] could 
provide only rough estimates of E^^min and -BpsF- Additional uncertainty is introduced in the 
estimate of -Bpsp by the uncertainty of the measurements of EBL which result in the uncertainty 
of the optical depth r in Eq. 1 . Further uncertainty is introduced in the derivation of E^^min from 
non-simultaneous data in GeV and TeV bands. Both HESS and Fermi measurements refer to the 
source spectra averaged over year(s) time scale. The reported HESS observations of the sources 
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took place in the period from 2005 to 2006 (22, 23, 36), whereas Fermi measurements were 
taken in 2008 and 2009. Up to now, no long-term variability was found in HESS observations 
(22, 23, 36, 43). 

The mean free path of primary multi-TeV gamma rays is of the order of ~ 80k Mpc. The 
largest structures in the universe, galaxy clusters, have typical sizes of the order of several Mpc 
and their volume filling factor is small. Most of the volume of the sphere of radius ~ 80k Mpc 
around lES 1 101-232, lES 0229-1-200 and lES 0347-121 is occupied by the voids in the Large- 
Scale Structure. 

Evidence for existence of magnetic fields in the voids provides a strong argument in favor 
of a cosmological origin of the fields serving as seeds for subsequent amplification in galaxies 
and galaxy clusters. Weak magnetic fields produced in the early universe are expected to fill 
the whole universe, including the voids. Contrarily, in the astrophysical models, the weak seed 
fields are created locally in (proto)galaxies and the field outside these structures should be close 
to zero. 

Cosmological magnetogenesis models consider generation of magnetic fields with a corre- 
lation length that does not exceed the size of the cosmological horizon and with energy density 
that does not exceed the critical density of the universe at the moment of magnetogenesis. Four 
broad classes of cosmological magnetogenesis scenaria are considered: magnetogenesis at the 
epoch of inflation, at the electroweak phase transition, at the epoch of quantum chromodynam- 
ics (QCD) phase transition and at the epoch of recombination {1-3, 44-51). The lower bound 
reported here excludes substantial parts of allowed parameter space for all the classes of cos- 
mological magnetogenesis models (Fig. [2l). 
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log(X, [Mpc]) 

Fig. 2: Light, medium and dark grey: known observational bounds on the strength and correla- 
tion length of EGMF, summarized in the Ref. (25). The bound from Big Bang Nucleosynthesis 
marked "BBN" is from the Ref. (2). The black hatched region shows the lower bound on the 
EGMF derived in this paper. Orange hatched regions show the allowed ranges of B, \b for 
magnetic fields generated at the epoch of Inflation (horizontal hatching) the electroweak phase 
transition (dense vertical hatching), QCD phase transition (medium vertical hatching), epoch of 
recombination (rear vertical hatching) (25). White ellipses show the range of measured mag- 
netic field strengths and correlation lengths in galaxies and galaxy clusters. 
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Supporting online material 

The FermifLAI data were filtered using tlie gtselect tool. Spectral analysis was done with 
the help of the likelihood technique using the gtlike tool, as explained in the Fermi/LAT data 
analysis threads (52). We have included point sources visible in FermifLAI images, Galactic 
diffuse emission and isotropic diffuse 7-ray emission components to the 7-ray emission model 
for the likelihood analysis. The mapcube file gllJem_v02.fits was used for the Galactic diffuse 
emission modeling, together with a corresponding tabulated model for the isotropic diffuse 
emission. The likelihood analysis for each object was performed inside a circular region with 
an angular radius of 14 degrees, centered on the source. 

To derive an upper limit on the source fluxes from the LAT data, we applied the following 
procedure. Assuming that the source flux in the 0.1-100 GeV band is a powerlaw dN^/dE = 
Nq {E/0.1 GeV)^'" we found the upper bound on normalization factor, Nq, for each T. To 
do this, we calculated the Test-Statistics (TS) (53) value as a function of A^^q and T. For each 
fixed value of T we plotted the TS value as a function of flux normalization to derive the 95% 
confidence level upper bound A0.95 (r) on the normalization factor, following procedure of Ref. 
(53). The resulting upper limits on the source fluxes shown as grey curves in Fig. 1 are the 
envelopes of the entire set of spectra dN^/ dE^ = A^o,95(r) {E^/ 0.1 GeV)~'". 

In order to calculate the bounds on the parameters of initial source spectra, we have fitted 
HESS data points in the 0.1-10 TeV energy band with model spectra for different Ecut and F, 
calculated under the assumption of zero magnetic field along the line of sight. For each model 
spectrum with given Ecut, T we found the normalization by fitting the model to the HESS data. 
Next, for each pair Ecut, T we found the of the fit and the total 7-ray flux emitted at the 
energies E > 0.1 TeV. The dependence of the and of the total flux on Ecut, T for the model 
fits of the HESS data is shown in Fig. [3l 
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Fig. 3: Grayscale: 68 and 95% confidence levels for the cut-off energy Ecut and photon index 
r found from the fitting of HESS spectra. Crosses mark the best- fit parameter values. Contours 
show the levels of integral energy flux above 0.1 TeV for the fitted spectra at each Ecut and 
r, with the increments of 2.5 x 10~^^ erg/cm^s. Arrows show the 95% confidence level lower 
bounds on the cut-off energy for the photon index F > 1.5 (upper arrows) and absolute lower 
bounds on Ecut (lower arrows). 
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r > 1.5 


aiir 


Source 




Ecut,min E^^min Bpsf{^ ^ -Dg) 


lES 0229+200 
lES 0347-121 
lES 1101-232 
H 2356-309 


3.8 TeV 39 GeV 3 x IQ-i^ q 
0.8 TeV 6 GeV 2 x IQ-^^ G 
1.0 TeV 
0.3 TeV 


1.2 TeV 30 GeV 3 x IQ-i^ q 
0.5 TeV 6 GeV 2 x 10^^^ G 
0.5 TeV 
0.14 TeV 



Table 1: Parameters of the model spectra and limits on EGMF for the analyzed sources. 

To find the values of parameters of initial source spectrum which minimize the cascade 
contribution to the source flux at zero EGMF strength, we chose a pair of cut-off energy Ecut,min 
and powerlaw index F lying in the 95% confidence contour and corresponding to the minimal 
integral energy flux above 0. 1 TeV. We consider the values of F softer than F = 1.5, a restriction 
commonly adopted in the modeling of blazar spectra in the very high-energy 7-ray band, see 
e.g. Refs. (22-24). We have verified that the results of our analysis (the existence of a lower 
bound on the EGMF strength) do not change if we relax the constraint F > 1.5 (see Table SI). 
The lower bounds on the cut-off energy found from the analysis of HESS spectra are given in 
Table SI, together with the estimates of E^^rnin and -Bpsf- 
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